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SUjMMARY 


Fracture  mechanics  has  been  found  to  be  a  useful  tool  for  understanding 
composite  delamination.  Analyses  for  calculating  strain  energy  release  rates 
associated  with  delamination  growth  have  been  developed#  These  analyses 
successfully  characterized  delaroination  onset  and  growth  for  particular 
sources  of  delamination.  Low  velocity  impact  has  been  found  to  be  the  most 
severe  source  of  composite  delamination.  A  variety  of  test  methods  for  mea¬ 
suring  interlaminar  fracture  toughness  are  being  developed  to  identify  new 
composite  materials  with  enhanced  delamination  resistance. 

INTRODUCTION 

As  high  strength,  low  density  fiber  reinforced  polymer  matrix  composites 
are  considered  for  highly  strained  primary  aircraft  structures,  increased 
attention  is  being  devoted  to  the  understanding  and  characterization  of  sig¬ 
nificant  composite  failure  modes.  Perhaps  the  most  commonly  observed  failure 
mode  in  composites  is  delamination,  \diich  is  a  separation  of  the  individual 
plies.  Delamination  may  result  from  low  velocity  impact,  from  eccentricities 
in  structural  load  paths  that  induce  out-of-plane  loads,  or  from  discontinu¬ 
ities  in  the  structure  which  create  local  out-^f -plane  loads.  Because  delam¬ 
ination  occurs  so  frequently,  interest  in  analyzing  and  characterizing  delam¬ 
ination  behavior  has  increased.  Fracture  mechanics  has  been  found  to  be  a 
useful  tool  for  this  purpose.  Summaries  of  some  of  the  significant  recent 
studies  on  delamination  are  included  in  this  review. 

A  Perspective 

There  is  a  fundamental  difference  between  the  way  linear  elastic  fracture 
mechanics  is  being  applied  to  de lamination  in  composites  and  the  way  it  is 
used  for  metals.  In  metals,  the  largest  crack  in  a  component  is  assumed  to 
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control  strength.  Failure  of  the  structure  may  be  directly  related,  by  the 
material's  fracture  toughness,  to  catastrophic  propagation  of  this  crack. 
However,  composite  delamination  in  structures  subjected  to  in-plane  loading  is 
a  subcritical  failure  mode  whose  effect  may  be  (1)  a  stiffness  loss  that  is 
benign  in  terms  of  the  structural  failure,  (2)  a  local  tensile  strain  concen¬ 
tration  in  the  load  bearing  plies  that  causes  tensile  failure,  or  (3)  a  local 
instability  that  causes  further  growth  which  ends  in  compressive  fed. lure.  In 
the  latter  two  cases  the  delamination  leads  to  a  redistribution  of  structural 
load  paths  which,  in  turn,  precipitates  structural  failure.  Hence,  delamina¬ 
tion  is  indirectly  responsible  for  the  final  failure  of  a  composite,  whereas* 
crack  propagation  in  a  metal  is  directly  responsible  for  structural  failure. 
Consequently,  interlaminar  fracture  toughness  for  composite  delamination  does 
not  have  the  same  significance  to  design  as  fracture  toughness  does  for  a 
metal.  The  designer  must  be  able  to  determine  the  consequence  of  the  delami¬ 
nation  in  terms  of  structural  load  redistribution  and  relate  this  to  some 
appropriate  structural  failure  criterion.  Nevertheless,  fracture  mechanics  is 
a  useful  tool  for  understanding  the  mechanics  of  delamination,  to  determine 
the  variables  that  control  when  delaminations  form  and  grow,  and  for  charac¬ 
terizing  the  inherent  delamination  resistance  of  the  composite.  This  paper 
focus  on  these  aspects  of  the  fracture  mechanics  characterization  of 
delamination  in  composite  materials. 

MECHANICS  OF  DELAMINATION 

Strain  Energy  Release  Rate  Analysis 

Some  potential  delamination  sites  are  the  common  design  details  that 
result  in  discontinuities  in  the  load  path.  These  discontinuities  give  rise 
to  interlaminar  stresses  even  when  the  remote  loading  is  in— plane  [1 J .  Some 
of  these  typical  design  details  (fig.  1)  are  (1)  a  straight  edge  with  a 
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stress-free  boundary  (often  called  a  "free  edge"),  (2)  a  curved  edge  such  as  a 
hole  boundary,  (3)  a  drop  off  of  the  interior  plies  of  a  laminate  to  taper 
thickness,  (4)  a  bonded  or  co-cured  joint,  and  (5)  a  bolted  joint.  Of  these 
cases,  the  free— edge  delamination  problem  has  been  studied  most  extensively, 
and  illustrates  the  benefits  of  fracture  mechanics  for  characterizing 
delamination  behavior. 

Several  investigators  have  analyzed  the  strain  energy  release  rate,  G, 

associated  with  edge  delamination  growth  [2,3,4].  Prom  finite  element  analysis 

they  have  found  that  once  the  delamination  progresses  beyond  a  distance  equal 

to  a  few  ply  thicknesses  from  the  edge,  G  reachehi'a  constant  plateau  given  by 

the  first  equation  in  figure  2.  This  equation,  derived  from  a  rule  of  mixtures 

and  laminated  plate  theory,  shows  that  G  is  independent  of  delamination  size 

and  depends  only  on  the  remote  strain,  e,  the  laminate  thickness,  t,  and  two 

modulus  terms,  and  B*,  that  correspond  to  the  laminate  modulus  before 

and  after  delamination.  In  Reference  [4],  the  critical  strain,  e^,,  measured 

at  the  onset  of  edge  delamination  in  the  -30/90  interfaces  of  eleven-ply 

[± 30/± 30/90/90 laminates  was  substituted  into  the  edge  delamiiiation  equation 

in  figure  2  to  determine  a  critical  G^,  for  TRornel  300/Narmco  5208* 

(graphite/epoxy)  composites.  This  Gq  was  then  used  to  predict  onset  strains 

for  edge  delamination  in  the  0/90  interfaces  of  [+45  /-45  /O  /90  ]  laminates 

n  n  n  ns 

of  the  same  material,  where  n  »  1,2,3  corresponded  to  and  8-,  16-,  and  24-ply 
laminate.  Figure  3  shows  that  the  predictions  agreed  well  with  experimental 
data,  and  captured  the  trend  of  decreasing  delamination  onset  strain  with 
increasing  thickness.  An  analysis  based  on  critical  values  of  interlaminar 
stresses  would  not  account  for  this  thickness  dependence  because  the 
interlaminar  stresses  are  independent  of  n  [3],  No  rigorous  analytical 

explanation  has  been  presented  to  explain  the  ability  of  the  strain  energy 

*The‘  use  of  trade  names  in  this  article  does  not  consitute  endorsement, 
either  expressed  or  implied,  by  the  National  .Aeronautics  and  Space 
Administration  or  AVSCOM.  , 


release  rate  (which  is  rigorously  defined  only  when  a  crack  of  finite  size 

exists)  to  predict  delamination  onset  when  a  stress  criterion  cannot. 

Although  the  interlaminar  stresses  in  the  interfaces  that  delaminate  are 

singular  at  the  straight  edge,  these  singularities  are  not  the  same  as 

those  that  exist  when  the  delamination  has  formed  [5,6].  Perhaps  the  most 

logical  explanation  is  the  presumption  that  inherent  interlaminar  flaws  of 

microscopic  dimension  are  initially  present  in  the  interface  near  the  edge, 

or  they  are  created  by  the  singular  stresses  at  the  edge,  but  they  do  not 

form  delaminations  of  a  finite  detectable  size  until  a  critical  G  value 

is  reached  for  delamination  extension  [4,7,8]. 

... 

A  similar  analysis  for  a  local  delamination  growing  from  a  matrix  ply 
crack  in  an  off-axis  ply  yielded  the  second  equation  shown  in  figure  2 
[9].  This  equation  also  shows  that  G  is  independent  of  delamination  size 
and  depends  only  on  the  remote  load,  P,  the  laminate  width,  W,  the  lami¬ 
nate  thickness,  t,  the  laminate  modulus,  E^am,  and  the  thickness  and 
modulus  of  the  uncracked,  locally  delaminated  regions,  t^o  and  En,, 
respectively.  Figure  4  shows  predictions  (open  symbols)  of  delamination 
onset  in  the  -25/90  interfaces  of  [±25/90n]g  laminates.  Under  tensile 
loading,  these  laminates  undergo  a  transition  from  edge  delamination 
(n  <  3)  to  local  delamination  growing  from  matrix  cracks  (n  >  4) .  A  criti¬ 
cal  Gq  was  calculated  from  the  n  =  2  edge-delamination  data,  and  then 
was  substituted  into  the  appropriate  equations  to  predict  delamination  on¬ 
set  strains  for  the  other  layups.  As  shown  in  figure  4,  the  G  anlyses 
predict  the  onset  of  both  types  of  delamination  well.  Furthermore,  just  as 
the  dependence  for  edge  delamination  was  predicted  by  G  analysis  (fig. 

3),  so  the  trend  of  decreasing,  e^,  with  increasing  thickness  of  the 
cracked  ply,  and  hence  increasing  n,  was  predicted  for  local  delamination 
(fig.  4). 


4 


The  ability  of  the  strain  energy  release  rate  to  correlate  delamination 
behavior  from  several  sources  and  account  for  volumetric  effects  in  the  form 
of  thickness  dependence  is  the  primary  reason  for  adopting  fracture  mechanics 
to  analyze  delaminatlon  problems.  Historically,  this  same  capability  for 
predicting  crack  growth  in  large  metallic  structures  based  on  fracture  tough¬ 
ness  measurements  of  small  test  specimens  led  to  the  acceptance  and  growth  of 
fracture  mechanics  in  the  1950's  [10]. 

Mixed-Mode  Behavior 

Because  a  delamination  is  constrained  to  grow  between  individual  plies, 
both  interlaminar  tension  and  shear  stresses  are  conunonly  present  at  the 
delamination  front.  Therefore,  delamination  is  often  a  mixed-mode  fracture 
process.  A  boundary  value  problem  must  be  formulated  and  solved  to  determine 
the  Gj',  Gjj,  and  Gjjj  components  of  the  total  strain  energy  release  rate. 

A  virtual  crack  extension  technique  has  been  combined  with  displacement-based 
finite  element  analysis  to  calculate  the  various  G  components  [21.  This 
technique,  in  conjunction  with  quasi-three-dimensional  finite  element  anal¬ 
yses,  was  applied  to  the  free-edge  de lamination  problem  by  several  authors 

[2.3.4.7.11.12.13] .  Figure  5  shows  the  result^  of  such  an  analysis  for  the 
[±30/±30/90/90]g  laminate  when  each  of  the  -30/90  interfaces  are  assumed  to 
delaminate.  As  noted  earlier,  the  total  G  raises  to  a  constant  plateau  when 
the  delamination  progresses  beyond  a  few  ply  thicknesses  from  the  edge. 

Figure  5  indicates  that  the  total  G  is  equal  to  the  sum  of  Gj  and  G^j 
components.  The  Gjjj  component  calculated  was  negligible.  Furthermore, 
like  the  total  G,  the  Gj  and  Gjj  components  are  also  independent  of 
delamination  size  for  edge  delamination. 

Similar  analyses  have  been  performed  for  other  composite  layups 

[3.11.12.13] .  Figure  6  shows  the  Gj  percentage  of  the  total  G  for  three 
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quasi-isotropic  layups  containing  ±45®,  0“,  and  90®  plies.  Each  layup  had  90® 
plies  in  the  center,  but  the  order  of  the  +45®,  -45®,  and  0®  plies  was  permu- 
tated.  Delaminations  were  modeled  between  the  90®  plies  and  adjacent  angle 
plies.  The  total  G  was  nearly  identical  for  the  three  layups  but,  as  shown 
in  figure  6,  the  Gj  component  varied  significantly. 

Because  of  the  complex  mixed-mode  nature  of  composite  delamination,  no 
closed  form  solutions  have  been  developed  to  lay  a  strong  theoretical  founda¬ 
tion  for  understanding  the  governing  parameters  that  control  de lamination 
behavior.  Recently,  however,  a  singular  integral  equation  approach  based  on 
®uisotropic  composite  laminate  elasticity  theory  was  used  to  establish  some  * 
general  asymptotic  solutions  for  edge  delamination  [14].  A  partially  closed 
crack  model  for  de lamination  was  introduced  to  remove  the  restriction  of 
physically  inadmissible  oscillatory  stress  singularities.  Both  the  global 
response  in  terms  of  strain  energy  release  rates  (Gj,  G^j,  Gjjj,  and 
*^T0TAL^  asymptotic  field  solution  in  terms  of  crack-tip  stress  inten- 

sity  factors  (Kj,  ^xi'  ®nd  ^xil^  were  calculated.  Anisotropic  laminate 
elasticity  solutions  were  obtained  for  relatively  simple  cases  of  delamination 
in  (+6/-6)g  laminates  [14].  Finite  element  methods  based  on  a  combined 
laminate  elasticity  solution  and  singular  hybrid  finite  element  formulation 
were  employed  for  fiber  composites  with  complex  lamination  variables  [15]. 
Figure  6  shows  the  G^  percentage  of  the  total  G  for  edge  delamination  in 
the  (±45/0/90 )g  family  analyzed  previously.  The  same  trend  in  Gj/G  with 
layup  is  apparent,  indicating  that  the  simple  approximate  techniques,  used 
with  the  more  common  displacement- based  finite  elements,  may  be  reasonably 
accurate  for  predicting  delamination  behavior.  The  virtual  crack  extension 
method  has  been  successfully  used  with  a  3D  NASTRAN  finite  element  analysis  to 
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predict  delamination  growth  from  both  a  manufacturing  flaw  and  an  open  hole  in 
the  F-16  composite  horizontal  tail  [l]. 

The  fracture  mechanics  analyses  of  other  sources  of  delamination  will 
continue  to  be  a  fruitful  research  topic  as  composite  materials  become 
candidates  for  primary  aircraft  structure.  Both  the  U.S.  Air  Force  and  the 
O.S.  Army  are  developing  damage  tolerance  and  durability  criteria  for  future 
composite  structures  which  may  ultimately  require  such  analyses  for  structural 
certification. 

DELAMINATION  GROWTH 

Delaminations  in  composite  laminates  may  interact  with  other  damage 
mechanisms  that  result  in  different  growth  behavior  than  may  be  anticipated 
from  an  elastic  analysis.  For  example,  under  quasi— static  loading,  edge 
delamination  formation  is. followed  by  matrix  cracking  in  the  90®  plies.  These 
matrix  cracks  may  extend  throughout  the  laminate  width,  ahead  of  the  delamina¬ 
tion  front  [4],  causing  a  perturbation  in  the  three-dimensional  stress  field 
at  the  delamination  front.  This  interaction  between  matrix  cracks  and  delami¬ 
nation  causes  stable  delamination  growth.  A  two-dimensional  linear  elastic 
fracture  mechanics  analysis  cannot  model  this  ^'three-dimensional  damage,  but  it 
may  be  used  for  semi -empirical  characterization  of  stable  growth  using  the 
R-curve  method  [16].  Figure  7  shows  an  R-curve  for  stable  edge  delamination 
growth  in  a  [±45/0/90] g  T300/5 208  laminate  [17].  Unlike  the  G  predicted 
from  an  elastic  analysis  that  is  independent  of  delamination  size,  a,  the 
delamination  resistance,  Gj^,  increases  with  a  as  the  matrix  cracking 
accumulates  ahead  of  the  delamination  front.  As  shown  in  figure  8,  this 
R-curve  was  used  to  predict  stable  edge-de lamination  growth  in  a  thicker 

[+452/-452/02/902]g  laminate,  and  in  a  [±30/±30/90/90]  laminate  with  a  dif- 

s  * 

ferent  layup. 
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If  the  composite  has  a  tough  majtrix,  however,  the  matrix  cracking  may 
be  suppressed.  Delaminations  will  form  at  higher  G  values,  but  they  will 
immediately  extend  through  the  laminate  width  in  an  unstable  manner  [7]. 
This  edge  delamination  behavior  in  a  toughened-matrix  composite  is  the 
opposite  of  the  stable  crack  growth  observed  in  a  ductile  metal  where  large 
plastic  zones  at  the  crack  tip  result  in  more  stable  growth. 

Another  source  of  delamination  resistance  that  results  in  stable 
growth  is  the  fiber  bridging  that  occurs  when  a  delamination  is  forced  to 
grow  between  plies  of  similar  orientation.  This  often  occurs  in  the  double 
cantilever  beam  (DCB)  test  used  to  measure  Mode  I  interlaminar  fracture 
toughness  [18,19,20].  Figure  9  shows  the  DCB  specimen,  vdiich  typically 
consists  of  24  plies  of  unidirectional  material.  As  shown  on  figure  9,  the 
specimen  compliance  is  a  non-linear  function  of  delamination  size.  In  ref¬ 
erence  [20],  a  Teflon*. coated  glass  fabric  was  placed  at  the  midplane  in 
the  end  of  the  beam  to  simulate  a  crack-like  sharp  notch.  Measurements  of 
Gjc  were  taken  at  the  initial  insert  length  {delamination  length  =  0)  and 
along  the  length  of  the  beam  as  the  delamination  grew.  Figure  10  shows  the 
R-curve  generated  from  a  DCB  specimen  of  Hercules  HMS/3501-6*  (graphite/ 
epoxy).  The  three-fold  increase  in  Gjq  from  toe  initiation  at  the  crack 
starter  until  values  stabilized  after  the  delamination  grew  over  one  third 
of  the  beam  legnto  was  attributed  to  fibers  debonding  from  the  matrix, 
bridging  the  delamination  plane,  flexing,  and  breaking.  A  similar  R-curve 
generated  for  an  Hercules  AS-1/3501-6*  DCB  specimen  showed  a  higher  Gj^ 
at  initiation  and  less  increase  in  Gj^  as  the  delamination  grew  (fig. 

10).  In  this  case,  much  less  fiber  bridging  was  observed  due  to  the 
stronger  interfacial  bond  between  the  AS-t  fiber  and  the  3501-6  matrix. 

Unlike  the  stable  delamination  growth  that  occurs  under  guasi— static 
loading  at  G  levels  above  G^,,  cyclic  de lamination  growth  may  occur  at 
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levels  well  below  the  static  interlaminar  fracture  toughness.  Figure  11  shows 
results  for  cyclic  edge  delamination  tests  at  maximum  cyclic  strain  levels 
below  the  delamination  onset  strain  measured  in  the  static  test  [11].  Delami¬ 
nations  formed  at  these  lower  cyclic  strains  after  a  certain  number  of  cycles, 

N.  values,  calculated  from  maximum  cyclic  strains  and  plotted  as  a  func¬ 
tion  of  the  number  of  cycles  to  delamination  onset,  are  shown  in  figure  11. 

These  values  drop  sharply  in  the  first  200,000  cycles  until  they  reach  a 

plateau  tantamount  to  a  threshold  for  delamination  onset  in  fatigue.  An 
elastic  analysis  of  the  [±35/0/90] ^  and  [0/±35/90]g  laminates  yields  an  iden¬ 
tical  total  G  for  a  given  nominal  strain,  but  the^  Gj  components  are  very 
high  and  low,  respectively.  As  shown  in  figure  11,  the  G^  for  these  two 
layups  measured. in  static  tests  (N  =  0)  are  different.  However,  the  G^, 
thresholds  under  cyclic  loading  are  nearly  identical.  Hence,  the  total  ‘ G 
does  not  determine  the  delamination  onset  under  static  loading,  but  does 
appear  to  govern  the  de lamination  onset  in  fatigue.  Furthermore,  comparison 
of  the  relatively  tough  Hexcel  HX205  matrix  composites  to  the  brittle  5208  matrix 

composite  shows  a  significant  improvement  in  the  static  G^,  yet  the  magnitude 
of  this  improvement  for  the  G^,  threshold  in  fatigue  is  much  less. 

The  comparisons  shown- in  figure  11  are  useful  for  selection  of  more 
durable  composite  matrices.  However,  for  damage  tolerance  analysis,  a  com¬ 
plete  characterization  of  delamination  growth,  from  the  onset  threshold  under 
cyclic  loading,  to  the  static  interlaminar  fracture  toughness,  G^,  is 

required  (fig. -12).  Several  investigators  have  used  the  mode  I  DCB  and  mixed¬ 
mode  cracked- lap-shear  (CLS)  specimens  to  characterize  this  delamination 

growth  [1,21,22,23].  They  have  developed- de lamination  growth  laws  of  the  form 
da  b  . 

—  =  cG  similar  to  cyclic  crack  growth  laws  used  to  characterize  fatigue 
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crack  growth  in  metals.  These  baseline  data  were  used  with  finite  element 
analyses  to  predict  cyclic  delamination  growth  in  structural  components  [1]. 

One  source  of  composite  delamination  growth  is  the  local  instability  that 
arises  from  the  buckling  of  a  delaminated  region  under  compression  loading 
[24,25].  Figure  13  shows  a  comparison  of  measured  delamination  cyclic  growth 
rates,  da/dN,  and  strain  energy  release  rate  components,  Gj  and 
calculated  from  a  geometrically  nonlinear  finite  element  analysis  for  a 
compressively  loaded  laminate  containing  a  through-width  delamination  (251. 
Both  Gj  and  the  delamination  growth  rate  increase  and  then  decrease  rapidly 
with  delamination  extension.  However,  Gjj  increase^^monotonically  to  a 
constant  level  with  delamination  extension.  Hence,  it  appears  that  tSj  is 
primarily  responsible  for  driving  the  delamination  in  the  buckled  region,  even 
though  the  Gjj  component  was  an  order  of  magnitude  larger. 

Figure  13  also  shows  the  tendency  for  the  delamination  cyclic  growth  rate 
to  decrease  with  delamination  length.  Hence,  the  growth  of  a  single  delamina— 
tion  under  cyclic  compression  may  arrest  and  hence  may  not  severely  degrade 
performance.  Figure  14  shows  plots  of  cycles  to  failure  as  a  function  of 
maximum  cyclic  strain  for  a  38-ply  T300/5208  ( graph! te/epoxy)  laminate,  typi¬ 
cal  of  an  upper  surface  transport  wing  skin,  subjected  to  cyclic  compressive 
loading  [26]  where  (1)  the  laminates  contain  a  single  circular  insert,  of 
various  diameters,  located  at  different  interfaces  to  simulate  delamination, 

(2)  the  laminates  have  been  subjected  to  low  velocity  impacts  over  a  range  of 
impact  energies,  and  (3)  the  laminate  contains  a  1.27  cm  (0.5  in)  diameter 
hole.  The  impacted  laminates  suffered  severe  degradation  in  the  compressive 
strain  at  failure.  In. most  cases  the  static  compressive  failure  strain  after 
impact  was  less  than  the  residual  compressive  strain  for  the  implanted  delami¬ 
nation  after  10®  fatigue  cycles. 
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Figure  15  shows  a  plot  of  cycles  to  failure  as  a  function  of  stress 
amplitude  for  [0/90/0/±45/0]g  graphite  epoxy  laminates  subjected  to  fully 
reversed  cyclic  loading^  either  in  the  undamaged  state^  or  following  an  impact 
with  a  potential  energy  per  unit  thickness  of  1790  J/m  [27].  All  failures 
occurred  in  the  compression  cycle.  The  data  in  figure  15  indicate  that  the 
compressive  strength  after  impact  is  very  low  compared  to  the  fatigue  behavior 
of  the  virgin  laminate.  Furthermore,  most  of  the  strength  reduction  occurs 
the  impact,  with  very  little  degradation  due  to  subsequent  cyclic 
loading. 

The  results  in  figures  14  and  15  indicate  that,  unlike  metals  where 
cyclic  loading  poses  the  most  severe  threat  of  crack  growth  leading  to  fail¬ 
ure,  the  most  severe  source  of  delamination  in  graphite  epoxy  structural 
laminates  is  low  velocity  impact.  Such  impacts  may  create  significant  delami- 

i  . 

nation  throughout  the  thickness  of  the  laminate,  degrading  compressive 
strength*  Furthermore,  this  extensive  delamination  may  occur  on  the  laminate- 
interior  with  little  or  no  visible  damage  apparent  on  the  laminate  surface 
[26-30]*  Some  research  has  just  begun  to  relate  the  strain  energy  release 
rate  for  delamination  growth  due  to  out-of -plane  loading  during  impact  to  the 
growth  of  multiple  delaminations  that  form  throughout  the  laminate  thickness* 
One  analysis  is  based  on  large  deflection  plate  theory  for  circular  isotropic 
plates  containing  multiple  delaminations  [31].  The  influence  of  plate  thick¬ 
ness,  support  conditions,  and  matrix  toughness  on  the  onset  and  propagation  of 
de lamination  was  predicted*  Figure  16  shows  a  comparison  of  the  predicted  and 
measured  load-displacement  record  for  a  16-ply,  circular  quasi-isotropic 
plate,  25  mm  in  diameter,  subjected  to  a  transverse,  load*  Seven  symmetrically 
loaded  delaminations  were  assumed  to  grow  simultaneously*  Good  agreement  was 
found  for  the  initial  delamination  onset  at  a  load  of  1  kN,  and  for  the 
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load~dlsplacement  curve,  up  to  the  point  that  fiber  breakage  occurred  in 
the  back  surface  ply,  at  a  load  of  3  kN.  More  refined  analyses  that  model 
the  specific  damage  details  observed  during  low  velocity  impact  are  needed 
to  further  clarify  the  relationships  between  impact  loading  and  interlami¬ 
nar  fracture. 

CHARACTERIZATION  OP  DELAMINATION  RESISTANCE 

Because  of  this  severe  delamination  threat  posed  by  non-vis ible,  low 
velocity  impact,  there  has  been  increased  interest  in  developing  new  tough¬ 
ened  matrix  resins  for  graphite  reinforced  composites.  Several  investiga¬ 
tors  have  identified  a  correlation  between  improved  compressive  strength 
after  Impact  and  improved  interlaminar  fracture  toughness  for  composltses 
with  toughened  matrix  resins  [26,29,30].  These  observations  provided  moti¬ 
vation  for  developing  new  test  methods  for  measuring  interlaminar  fracture 
toughness. 

Many  investigators  have  used  a  double  cantilever*  beam  (DCB)  specimen 
(fig.  9)  to  measure  critical  values  of  the  crack  opening  ( inter launinar 
tension)  fracture  toughness,  Gjq  [18-23,32-34].  The  edge  delamination 
test  (EDT)  (fig.  17)  is  also  being  used  to  measure  the  mixed— mode  fracture 
toughness  of  composites  [7]  and  determine  fracture  mode  dependence  [35,- 
36].  Figure  18  shows  a  comparison  of  the  interlaminar  fracture  toughness 
of  graphite  composites  with  very  brittle  or  tough  matrix  resins  tested 
using  a  24-ply  unidirectional  DCB  as  well  as  three  EDT  layups  with  low,  in¬ 
termediate,  and  high  percentages  of  Gj.  The  solid  symbols  represent  mode 
I  component  of  measured  G^.  Results  indicated  that  only  the  Gj  compon¬ 
ent  controlled  interlaminar  fracture  for  the  brittle  matrix  composites 
(5208,  2202-1)*  whereas  both  interlaminar  tension  and  shear  contributed  to 
delamination  in  tough  matrix  composites  (HX205/P185) *.  Both  the  DCB  and 
EDT  tests  are  being  used  to  evaluate  new  composite  materials  developed  for 
heavily  loaded  large  composite  wing  skins  under  the  NASA  Aircraft  Energy 
Efficiency  (ACEE) ’ project  [37]. 
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Several  other  tests  are  also  being  used  to  measure  interlaminar  fracture 
toughness*  A  pure  mode  I  version  of  the  EDT  test  has  been  proposed  [38]#  The 
cracked- lap  shear- (CLS)  test  for  mixed  mode  has  been  used  by  several 

investigators  [19,21,22] •  Currently  an  ASTM  task  group  is  conducting  a  round 
robin  test  program  using  the  DCB,  EDT,  and  CLS  specimen  as  a  prelude  to  devel¬ 
oping  ASTM  standards.  Other  investigators  have  used  unidirectional  beams 

loaded  in  three-point  bending  to  estimate  g  [19,39,40,41]  or  loaded 

lie 

asymmetrically  to  measure  mixed  mode  values  with  different  ratios  of 

Gi/Gii  [19,40]. 

The  recognition  of  the  importance  of  interlaminar  fracture  toughness  to 
composite  damage  tolerance  has  led  to  the  development  of  new  toughened  matrix 
composites.  Because  resin  manufacturers  would  prefer  to  measure  the  fracture 
‘toughness  of  a.  small  neat  resin  specimen,  there  has  been  renewed  interest  in 
composite  micromechanics  analyses.  TTie  ultimate  goal  of  such  analyses  would 
be  to  relate  the  fracture  behavior  of  the  neat  resin  to  the  interlaminar  frac¬ 
ture  behavior  of  the  composite.  Such  analyses^would  guide  the  development  of 
neat  resin  fracture  tests  that  would  yield  realistic  estimates  of  improvements 
in  interlaminar  fracture  toughness  of  composites. 

CX3NCLUSIONS  AND  FUTURE  DIRECTIONS 

Although  the  promise  of  composite  materials  has  been  known  for  several 
decades,  the  application  of  these  materials  to  primary  aircraft  structure  is 
just  beginning.  As  suitable  toughened-resin  composites  are  developed  to 
overcome  the  existing  impact  problems,  methods  for  characterizing  their 
behavior  will  have  to  be  re-evaluated.  For  example,  it  may  become  necessary 
to  use  a  J-integral  analysis  to  correlate  delamination  behavior  between 
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different  layups  and  structural  discontinuities*  Furthermore,  fatigue  may 
reappear  as  the  limiting  factor  for  damage  tolerant  design.  In  any  event,  the 
ability  to  anticipate  and  quantify  delamination  failure  using  fracture  mechan¬ 
ics  analyses  will  still  play  an  important  role  in  the  understanding  of 
composite  tehavior,  in  establishing  the  damage  tolerance  and  durability  design 
criteria  for  certification,  and  in  the  development  of  improved  composite 


materials 
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Fig  2  Sources  of  delamination  In  unnotched  la*mlnates 
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Fig  3  Edg 
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Fig  4  Prediction 
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Fig  5  Strain  energy  release  rate  components  for  edge  delamination 
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Fig  6  Mode  I  percentage  of  total  G  for  edge  delamination  in  0/90  interfaces  of 
quasi-lsotropic  laminates 
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Fig  10  R -curves  for  delamination  growth  in  unidirectional  double  cantilever  beam  specimens 
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Fig  11  Interlaminar  G^,  as  a  function  of  fatigue  cycles 
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Fig  14  Cyclic  compression  fatigue  life  trends  for  T300/5208  graphite  epoxy  laminates 
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fully  reversed  cyclic  loading 
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Fig  17  Edge  delamination  test  for  mixed-mode  interlaminar  fracture  toughness 
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